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ABSTRACT: In this paper, we introduce an approach to increase the density of field-effect transistors 
framework a compact transimpedance amplifier. In the present theoretical framework, the approach we 
consider in the manufacturing the inverter in heterostructure with specific configuration. Several required 
areas of the heterostructure should be doped by diffusion (alternatively they could be doped by ion 
implantation). After the considered doping dopant and radiation defects should by annealed framework 
optimized scheme. In this study, we also consider an approach to decrease the values of mismatch-induced 
stress in the considered heterostructure. An analytical approach is used to analyze the mass and heat 
transport in heterostructure during the manufacturing of integrated circuits with account of mismatch- 


induced stress. 
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induced stress and porosity of materials. 


https://doi.org/10.29294/IASE.8.3.2022.2242-2251 


INTRODUCTION 

In the present days, several problems of the 
solid state electronics (such as increasing of 
performance, reliability and density of elements of 
integrated circuits: diodes, field-effect and bipolar 
transistors) are intensively solving [1-6]. To 
increase the performance of the devices, it is 
attracted an interest in the determination of 
materials with higher values of charge carriers 
mobility [7-10]. One way to decrease the 
dimensions of elements of integrated circuits is by 
manufacturing them in thin film heterostructures 
[3-5,11]. In this case, it is possible to use 
inhomogeneity of heterostructure and necessary 
optimization of doping of electronic materials [12] 
and to develop an epitaxial technology to improve 
these materials (including the analysis of 
mismatch induced stress) [13-15]. An alternative 
approaches to increase the dimensions of 
integrated circuits are using of laser and 
microwave types of annealing [16-18]. 


Framework, in this paper we introduce an 
approach to manufacture field-effect transistors. 
This approach gives a possibility to decrease their 
dimensions with increasing their density 
framework a compact transimpedance amplifier. 
We also consider the possibility to decrease 
mismatch-induced stress to decrease quantity of 
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defects, generated due to the stress. Also, we 
consider a heterostructure which consist of a 
substrate and an epitaxial layer (see Fig. 1), anda 
buffer layer between the substrate and the 
epitaxial layer. The epitaxial layer includes itself 
several sections, which were manufactured by 
using another materials. These sections have been 
doped by diffusion or ion implantation to 
manufacture the required types of conductivity (p 
or n). These areas became sources, drains and 
gates (see Fig. 1). After doping, it is required to 
annealing of dopant and/or radiation defects. The 
objective of the present paper is to analysis the 
redistribution of dopant and radiation defects to 
determine under conditions, which correspond to 
decreasing of elements of the considered voltage 
reference and at the same time to increase their 
density. At the same time, we consider a 
possibility to decrease mismatch-induced stress. 


Substrate 


Fig. 1b Heterostructure with a substrate, epitaxial 
layers and buffer layer (view from side) 
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Fig. 1a. Structure of the considered amplifier [19] 


METHOD OF SOLUTION 


To solve our aim, we determine and analyzed the spatio-temporal distribution of concentration of dopant 
in the considered heterostructure. We determine the distribution by solving the second Fick's law in the 
following form [1,20-23] 


AC(x, y,z,t) 6 | ,OC(xmy,zt)|, 2 AC(x,y,z.t)| @| AC(x,y,z,t) 


= —|D +—|D 
Ot Ox Ox Oy Oy OZ OZ 
0 | Dy rT, 
+Q—| ZV u(x, y,z,t) (C(x, y,W.t)dW |+ 

Ox| kT 0 

0 | D, 7 
+Q—| SV u(x, y, zt) (C(x, y,W.t)dW | + (1) 

Oy| kT 0 


ai O Des 6 w(x, y,z,t) +4 O Des 6 w(x, y,z,t) 7" O Des 6 p(x, y,z,t) 
Ox|VkT Ox Oy|VkT Oy Oz|VkT Cz 


with boundary and initial conditions 


aC(x, y, z,t) 0C(x, y, zt) AC(x, y,z,t) 
—— = 0, ————+ =), —— 0-0 (x,V,2,0)=fc (xy,Z), 
Ox Hag Ox a Oy , 
AC (x, y,z,t) _0 AC (x, y, z,t) _¢ AC(x, y, z,t) _0 
0 x x=L, é < z=0 é < x=L, 
Here C(x,y,z,t) is the spatio-temporal distribution of dopant; Vs is the symbol of surficial gradient; 
of concentration of dopant; Q is the atomic volume 
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L, 
{C(x, y, z,t)d Z is the surficial concentration 
0 


of dopant on interface between layers of 
heterostructure (in this situation we assume, that 
Z-axis is perpendicular to interface between layers 
of heterostructure); s(x,y,z,t) and so(x,y,z,t) are 
the chemical potential due to the presence of 
mismatch-induced stress and porosity of material; 


De =D x y,2,T) Lg 


Ds Dye y, 2,7) 1+¢; 


Here Dz (x,y,z,T) and Dis (x,y,z,T) are the spatial 
(due to accounting all layers of heterostruicture) 
and temperature (due to Arrhenius law) 
dependences of dopant diffusion coefficients; T is 
the temperature of annealing; P (x,,z,T) is the 
limit of solubility of dopant; parameter y depends 
on properties of materials and could be integer in 
the following interval ye[1,3] [24]; V(x,y,z,t) is the 
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D and Ds are the coefficients of volumetric and 
surficial diffusions. Values of dopant diffusions 
coefficients depends on properties of materials of 
heterostructure, speed of heating and cooling of 
materials during annealing and spatio-temporal 
distribution of concentration of dopant. 
Dependences of dopant diffusions coefficients on 
parameters could be approximated by the 
following relations [24-26] 


V (x, y, zt V(x, y, zt 
l+¢, Side Mee ee ) 
V (x, y, Z,¢) V’ (x, y, Zt) 


I+¢, a +S (Cai - (2) 


spatio-temporal distribution of concentration of 
radiation vacancies; V is the equilibrium 
distribution of vacancies. Concentrational 
dependence of dopant diffusion coefficient has 
been described in details in [24]. Spatio-temporal 
distributions of concentration of point radiation 
defects have been determined by solving the 
following system of equations [20-23,25,26] 
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Ox|VkT Ox 


ae = (3) 
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Here p =I,V; I (x,y,z,t) is the spatio-temporal 
distribution of concentration of radiation 
interstitials; [* is the equilibrium distribution of 
interstitials; Di(x,y,z,T), Dv(x,y,z,T), Dis(xy, z,T), 
Dys(x,y,Z,T) are the coefficients of volumetric and 
surficial diffusions of interstitials and vacancies, 
respectively; terms V2(x,y,z,t) and [(x,y,z,t) 
correspond to generation of divacancies and 
diinterstitials, respectively (see, for example, [26] 
and appropriate references in this book); 
kiv(xy,Z,T), kiu(xy,z,T) and kyvv(x,y,z,T) are the 
parameters of recombination of point radiation 
defects and generation of their complexes; k is the 
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Boltzmann constant; w = a3, a is the interatomic 


distance; ¢ is the specific surface energy. To 
account porosity of buffer layers we assume, that 
porous are approximately cylindrical with average 


values Fr = Pe + y; and z1 before annealing 


[23]. With time small pores combine into large 
ones [27]. With time large pores became larger 
due to absorbing the vacancies and became more 
spherical [27]. Distribution of concentration of 
vacancies in heterostructure, existing due to 
porosity, could be determined by summing on all 
pores, i.e. 


VeoeneS Ss V(xtia,y+jB.ztky.t), Ra\ery ee. 


i=0j=0k=0 


Here a, Zand yare the average distances between centers of pores in directions x, y and z; |, m and nare the 


quantity of pores inappropriate directions. 


Spatio-temporal distributions of divacancies @y (x,y,z,t) and diinterstitials ® (x,y,z, t) could be determined 


by solving the following system of equations [25,26] 
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+k,(x, yz.) p(x, y,zt) (5) 
with boundary and initial conditions 
O® ,(x,y,z,t) = O® , (x, y,z,t) 0 O® , (x, y,z,t) =, 
Ox ba Ox _ Oy i _— 
O® ,(x, y,z,1) _0. O® ,(x,y,z,t) o O® , (x, y,z,t) i a 
O ¥ y=L, 0 < z=0 C2 z=L, 


Pp (x,y,2,0)=fap (x,y,2). 
Here Dai(x,y,z,T), Dov(x,y,Z,T), Das (x,y,z,T) and Devs(x,y,z,T) are the coefficients of volumetric and surficial 
diffusions of complexes of radiation defects; ki(x,y,z,T) and kv(x,y,z,T) are the parameters of decay of 


complexes of radiation defects. 


Chemical potential sm in Eq.(1) could be determine by the following relation [20] 
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M=E(Z)Qoy [ui(Xy,z,) +ui(xy,Z,0)]/2, (7) 


where E(z) is the Young modulus, oj is the stress tensor; Uj; = 


1 Ou, , Ou, 
2\ Ox; OX, 


is the deformation 


tensor; ui, uj are the components u:(x,y,z,t), uy(xy,z,t) and u:-(x,y,z,t) of the displacement vector 


u (x, y,Z, t); Xi, Xj are the coordinate x, y, z. The Eq. (3) could be transform to the following form 


u(x y,z,t)= 


o(z)6; Ou,(x, y, z,t) 
*  1-20(z) Cx 


where oa is Poisson coefficient; & = (as-azz)/az is 
the mismatch parameter; ds, az, are lattice 
distances of the substrate and the epitaxial layer; 
K is the modulus of uniform compression; / is the 
coefficient of thermal expansion; T; is the 


6u,(x, y,Z,t) | du, (x, y, z,t) 1 6u,(x, y, Zt) | du,(x, y, Zt) 
Ox, Ox, 


2 Ox, Ox, 


36; ~K(z)B(z)[T (x y.z,.t)-T | 6 = E(2) 


equilibrium temperature, which coincide (for our 
case) with room temperature. Components of 
displacement vector could be obtained by solution 
of the following equations [21] 


2 
(22 u(x, y, Zt) Z O0,,(%,¥.%t) | OO i Je8) : do,,(x, y,z,t) 


r ar ax ay Az 
2 
P ee Hj), z,t) Z do,,(x, y,Z,t) rn 00, (x, y, z,t) : Tom (x, y, Zt) s 
Ot Ox Oy OZ 
. 0 Veet 
P ie Ue (a yz) _ do.,(x, Y.Zst) o.,(x y,Z ). do_(x, y,z,t) 
Ot Ox Oy Oz 
where 
o, = E(z) du(x.y,z.t)  Oul%yzt) 4) du,(x, y, zt) + K(z)6, x 
' 2fl+o(z) Ox, Ox, a Om / 
x Ou, y zt) _ B(z)K (z)|T (x, y, 2t)—T.|, p (z) is the density of materials of 
X 
beeouace oy Is the Kronecker symbol. The conditions for the system in Eqs. (8) could be written in 
the form 
0 H(0,y,z,t)_ 9 O8(L,.¥Zt)_ 9) FH(40,z,t)__ 6 W(xL,.2t) 9 
Ox , Ox Oy Oy — 
d(x y.O.t) _ 9, g fi(xy Lt) _ 9 ii(x, y,z,0) =i, ; ii(x, y, 2,0) = thy. 
Oz OZ 


We calculate the spatio-temporal distributions 
of concentrations of dopant and radiation defects 
by solving the Eqs.(1), (3), (5) and (8) framework 
standard method of averaging of function 
corrections [28]. In this paper, we also determine 
the concentrations of dopant, concentrations of 


radiation defects and components of displacement 
vector by using the second-order approximation 
framework method of averaging of function 
corrections. This approximation is usually enough 
good approximation to make qualitative analysis 
and to obtain some quantitative results. All 
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obtained results have been 


DISCUSSION 


checked by 


In this section, we analyzed the dynamics of 
redistributions of dopant and radiation defects 
during annealing, under the influence of 
mismatch-induced stress and modification of 
porosity. Typical distributions of concentrations of 
dopant in heterostructures are presented on Figs. 
2 and 3 for diffusion and ion types of doping, 
respectively. These distributions have been 
calculated for the case, when the value of dopant 
diffusion coefficient in doped area is larger than in 
nearest areas. The Figures shown that the 
inhomogeneity of heterostructure gives us 
possibility to increase the compactness of 
concentrations of dopants and at the same time to 
increase homogeneity of dopant distribution in 
doped part of epitaxial layer. However, this 
approach of manufacturing of bipolar transistor it 


1 L,LyL. 
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comparison with results of numerical simulations. 


is necessary to optimize annealing of dopant 
and/or radiation defects. Reason of _ this 
optimization is following. If annealing time is 
small, the dopant did not achieve any interfaces 
between materials of heterostructure. In this 
situation one cannot find any modifications of 
distribution of concentration of dopant. If 
annealing time is large, distribution of 
concentration of dopant is too homogenous. We 
optimize annealing time framework recently 
introduces approach [29-37]. Framework this 
criterion we approximate real distribution of 
concentration of dopant by step-wise function (see 
Figs. 4 and 5). Farther we determine optimal 
values of annealing time by minimization of the 
following mean-squared error 


U=—— J { f[c(x.y.z.0)-w(x.y.2))dzd ydxx, (15) 
000 


LLL, 


C(x.0) 


Fig.2. Distributions of concentration of infused dopant in heterostructure from Fig. 1 in direction, which is 
perpendicular to interface between epitaxial layer substrate. Increasing of number of curve corresponds to 
increasing of difference between values of dopant diffusion coefficient in layers of heterostructure under 
condition, when value of dopant diffusion coefficient in epitaxial layer is larger, than value of dopant 


diffusion coefficient in substrate. 
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Fig.3. Distributions of concentration of implanted dopant in heterostructure from Fig. 1 in direction, which 
is perpendicular to interface between epitaxial layer substrate. Curves 1 and 3 corresponds to annealing 
time © = 0.0048(Lx2+L,2+Lz7) /Do. Curves 2 and 4 corresponds to annealing time © = 0.0057 (Lx2+Ly2+Lz?) /Do. 
Curves 1 and 2 corresponds to homogenous sample. Curves 3 and 4 corresponds to heterostructure under 
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condition, when value of dopant diffusion coefficient in epitaxial layer is larger, than value of dopant 


diffusion coefficient in substrate 


C(x,0) 


0 


L 


x 


Fig. 4. Spatial distributions of dopant in heterostructure after dopant infusion. Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing time. 
Increasing of number of curve corresponds to increasing of annealing time 


where yw (x,y,z) is the approximation function. 
Dependences of optimal values of annealing time 
on parameters are presented on Figs. 6 and 7 for 
diffusion and ion types of doping, respectively. It 
should be noted, that it is necessary to anneal the 
radiation defects after ion implantation. One could 
find the spreading of concentration of distribution 
of dopant during annealing. In the ideal case 


C(x,0) 


distribution of dopant achieves appropriate 
interfaces between materials of heterostructure 
during annealing of radiation defects. If dopant did 
not achieve any interfaces during annealing of 
radiation defects, it is practicably to additionally 
anneal the dopant. In this situation optimal value 
of additional annealing time of implanted dopant 
is smaller, than annealing time of infused dopant. 


0 


L 
x 


Fig. 5. Spatial distributions of dopant in heterostructure after ion implantation. Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing time. 
Increasing of number of curve corresponds to increasing of annealing time 
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Fig.6. Dependences of dimensionless optimal annealing time for doping by diffusion, which have been 
obtained by minimization of mean-squared error, on several parameters. Curve 1 is the dependence of 
dimensionless optimal annealing time on the relation a/L and €= y= 0 for equal to each other values of 
dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of dimensionless 
optimal annealing time on value of parameter ¢ for a/L=1/2 and €= y= 0. Curve 3 is the dependence of 
dimensionless optimal annealing time on value of parameter € for a/L=1/2 and ¢ = y= 0. Curve 4 is the 
dependence of dimensionless optimal annealing time on value of parameter yfor a/L=1/2 and e= €=0 


0.0 0.1 0.2 0.3 0.4 0.5 
a/L, &, € ¥ 


Fig.7. Dependences of dimensionless optimal annealing time for doping by ion implantation, which have 
been obtained by minimization of mean-squared error, on several parameters. Curve 1 is the dependence of 
dimensionless optimal annealing time on the relation a/L and €= y= 0 for equal to each other values of 
dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of dimensionless 
optimal annealing time on value of parameter ¢ for a/L=1/2 and €= y= 0. Curve 3 is the dependence of 
dimensionless optimal annealing time on value of parameter ¢ for a/L=1/2 and ¢ = y= 0. Curve 4 is the 
dependence of dimensionless optimal annealing time on value of parameter yfor a/L=1/2 and e= €=0 


1.0 


0.8 
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Fig. 8. Normalized dependences of component uz of displacement vector on coordinate z for nonporous 
(curve 1) and porous (curve 2) epitaxial layers 
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Fig. 9. Normalized dependences of vacancy concentrations on coordinate z in unstressed (curve 1) and 
stressed (curve 2) epitaxial layers 
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Further, we analyzed influence of relaxation of 
mechanical stress on distribution of dopant in 
doped areas of heterostructure. Under following 
condition «< 0 one can find compression of 
distribution of concentration of dopant near 
interface between materials of heterostructure. 
Contrary (at s>0) one can find spreading of 
distribution of concentration of dopant in this 
area. This changing of distribution of 
concentration of dopant could be at least partially 
compensated by using laser annealing [37]. This 
type of annealing gives us possibility to accelerate 
diffusion of dopant and another process in 
annealed area due to inhomogenous distribution 
of temperature and Arrhenius law. Accounting 
relaxation of mismatch-induced stress _ in 
heterostructure could leads to changing of optimal 
values of annealing time. At the same time 
modification of porosity gives us possibility to 
decrease value of mechanical stress. On the one 
hand mismatch-induced stress could be used to 
increase density of elements of integrated circuits. 
On the other hand could leads to generation 
dislocations of the discrepancy. Figs. 8 and 9 show 
distributions of concentration of vacancies in 
porous materials and component of displacement 
vector, which is perpendicular to interface 
between layers of heterostructure. 


CONCLUSION 


In this paper, we modeled the redistribution of 
infused and implanted dopants with account 
relaxation mismatch-induced_ stress during 
manufacturing field-effect § heterotransistors 
framework a compact transimpedance amplifier. 
We also formulated the recommendations for 
optimization of annealing to decrease dimensions 
of transistors and to increase their density and the 
recommendations to decrease mismatch-induced 
stress. Analytical approach to model diffusion and 
ion types of doping with account concurrent 
changing of parameters in space and time has 
been introduced. At the same time the approach 
gives us possibility to take into account 
nonlinearity of considered processes. 
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